OSU-95-101 



The Gluon Condensate and Running Coupling of QCD 

Hua-Bin Tang and R. J. Furnstahl 

Department of Physics 
The Ohio State University, Columbus, Ohio ^3210 
(February 16, 1995) 

Abstract 

An expression for the photon condensate in quantum electrodynamics is pre- 
sented and generalized to deduce a simple relation between the gluon conden- 
sate and the running coupling constant of quantum chromodynamics (QCD). 
Ambiguities in defining the condensates are discussed. The values of the gluon 
condensate from some Ansatze for the running coupling in the literature are 
compared with the value determined from QCD sum rules. 
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The running coupling constant of quantum chromodynamics (QCD) in the low- 
momentum region is so far scarcely known, and yet it is associated with various approaches 
to low-energy nonperturbative QCD. For example, in the Dyson-Schwinger equation (DSE) 
approach, one usually assumes some Ansatz for the running coupling constant in momen- 
tum space to approximate the gluon propagator. In quark potential models, one also needs 
to know the associated interaction potential. Thus, apart from the well-known constraints 
from confinement and asymptotic freedom, additional constraints on the running coupling 
constant in the low-momentum region would be very useful. One such constraint is discussed 
here. 

The phenomenological success of QCD sum rules || has stimulated many investigations 
|3|-U] on the validity of the operator-product expansion (OPE) |§. An important result 
has been the realization that one should introduce a renormalization scale [i in the Wilson 
coefficients as well as the condensates PJlO|. The high-momentum (p > \l with p being 
Euclidean) parts of the operator product are put into the coefficients and the low-momentum 
(p < yu) parts are parametrized by the condensates. Although the formal validity of the OPE 



is now generally accepted, there are still many discussions [pTT] |13]] concerning the practical 
separation of the low- and high-momentum contributions and the ambiguities associated 
with the infrared renormalon singularities |]Tl|,|T3|,|in| . 



A related issue is whether one can separate the perturbative contributions in the form of 
powers of [i from the so-called "genuine" nonperturbative effects in the condensates. While 
Ref. || suggests that the separation is possible, Ref. J?| seems to imply the opposite. In 
this Letter, we hope to clarify the situation by using some simple examples. We proceed to 
the discussion of the photon condensate in quantum electrodynamics (QED). Generalizing 
this QED result, we deduce a simple relation between the gluon condensate and the running 
coupling constant of QCD. Some Ansatze for the running coupling in the literature are then 
examined. 

We start with the simplest example of the composite operator cf 2 with <p(x) being a free 
scalar field. One can define <fi 2 as the x — > limit of the regular part of T[<p(x)<p(0)]: 



i> 2 = lim 



T[<f>(x)cf>(0)}-S(x) , (1) 



where T stands for the usual time-ordering and S(x) is the part of T[<p(x)<f)(0)] that is 
singular at x = and it can contain other composite operators with singular coefficients 
when interactions are present. In general S(x) cannot be uniquely defined since, at the very 
least, it may contain a regular function multiplied by the singular logarithm ln(/i 2 x 2 ), where 
fi is some arbitrary scale to make the argument dimensionless. We call S(x) minimal when 
the regular parts of S(x) come only from their multiplication with the singular function 
ln(/i 2 x 2 ) but not with powers of \x. It is instructive to compare the vacuum matrix element 
of <fi 2 in cutoff and dimensional regularizations with S(x) being minimal. 
Let's first consider using cutoff regularization (CR). We can write 



+ 



<l 1 [l + ik-x - ^(A;-x) 2 H ] . (2) 



fc<M (2vr) 4 k 2 + M 2 
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where M is the scalar mass and we have introduced the cutoff scale fi with // 3> M. The 
momentum is Euclidean. Notice that the first integral is not a minimal singular part. Indeed 
performing the integral term-by-term one sees that it contains the regular piece 

M 2 , u 2 ^ , a* u 2 \ o 

where 7 = 0.577215 is Euler's constant. The above smooth piece should be removed from 
the first integral to make it minimal. After such a removal and using Eq. ([I]), we obtain 

d 4 fc 1 M 2 - 2 

(2vr) 4 W- 
r2 ,,2 



<* 2|0> = Lwf - Obfe + 1 - 2 -> + ln4 ) 



^_( ln ^ + 1 _ 27 + lD4 ). (4) 

If we consider dimensional regularization (DR) instead, we can also make an analogous 
separation, 

4 _ » r d^ (_1)«M 2 " 



(O|T[0(x)0(O)]|O) = / . 4 - d E/7^p e 



n=0 ' 

00 



(2vr) d (jfc 



2^n+l 



^ ^ n!7 (2vr) rf k 2 + M 2 ' K ' 

where the first sum contains the singular piece JTB] . Eq (||) is justified as follows. Subtracting 
the first (singular) sum on the right-hand side of Eq (0) from the propagator on the left-hand 
side, expanding the exponential e lk ' x in the resulting regular expression, and using 



d d k{k 2 Y = (6) 

with arbitrary (3, one arrives at the second sum.[] The terms in the second sum of Eq ([|) are 
the so-called non-singular parts in the OPE and are proportional to (^-functions and their 
derivatives in momentum space . Notice again that the first sum in Eq (|5|) is not minimal. 
Performing the integrals and letting d — > 4, one finds that the following regular Taylor series 
(singular in d) 

M 2 ( 1 \ M 4 / 2 i 1 x 2 

=■(- 7-Iiitt 7 + l-ln7rx 2 H (7) 

167r 2 Vd-4 ' I 128tt 2 M-4 ' ' v ' 

should be removed to make the singular integral minimal. Eq. ([!]) then leads to 

4 H f d d k 1 M 2 / 2 , x 

M 2 11 2 



x It may be a helpful exercise for skeptical readers to explicitly check Eq (0) for d=l when all 
the integrals are trivial. In that case, the separation in Eq (0) is into a purely singular piece: 
-(M/2) En=o( M2x2 ) n /( 2n + !) ! and a regular piece: (1/2M) E^o^ 2 ^ 2 )™ /( 2ra ) ! - 
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which is the same as Eq. Thus the result for (O|0 2 |O) is independent of the regularization 
schemes. 

In relation to the usual definition of composite operators, the subtractions in Eqs. (|J) 
and (H) correspond to additive renormalization [I7l,j8| or mixing with the unit operator. The 



above simple exercise suggests that the mixing of composite operators with the poles in 
DR, or power terms (~ fi 4 , fi 2 ) in CR, are artifacts of defining the singular part S(x) to 
include these same terms. To further assure this point we note the following. The use of 
Eq. (|^) with DR in evaluating Feynman diagrams can be justified in some cases by a careful 
consideration of jacobian factors or extra counterterms proportional to a 5-function at the 
origin In the OPE with DR, the validity of Eq. (|6|) is further related to the fact that 
we apply the same regularization and neglect similar integrals in the Wilson coefficients. In 
the example of Eq. (|5|), the use of DR to regularize the infrared singularities in the first sum 
justifies the use of Eq. (Q) in calculating the integrals in the second sum. In other words, 
by using DR, we have implicitly and simultaneously neglected power divergent terms in the 
Wilson coefficients and the condensates. Thus neglecting the power terms in CR amounts 
to doing the same thing implicit in DR. Note that the power terms in the coefficients come 
from the infrared side of the momentum region, while they appear in the condensate from 
the ultraviolet side. 

In an asymptotically free theory such as QCD, the genuine infrared nonperturbative con- 
tributions, in contrast to the power terms, appear only in the condensates if the coefficient 
functions are calculated in perturbation theory as in the QCD sum-rule approach. Thus the 
nonvanishing gluon condensate is physically significant and not just a convenient definition, 



in contrast to Ref. |T^j. In the practical use of the OPE in QCD sum rules, Wilson coeffi- 
cients are calculated without the power dependence 0. As a result the phenomenological 
condensates as determined from the sum rules should not contain the power terms. In fact, 
as will be shown below, if the power terms were not excluded, they would be comparable to 
the phenomenological gluon condensate at fi ~ 1 GeV, in contradiction with the sum-rule 
assumption of the dominance of the genuine nonperturbative effects. Thus caution should 
be exercised when one uses such phenomenological condensates or makes comparison with 
them. 

One might suspect that the above free-field exercise is too simple to generalize, especially 
when there are infrared renormalon ambiguities. In fact, there seems to be disagreement 
in the literature ||[7| on whether one can separate the power terms from the genuine non- 
perturbative effects in the two-dimensional O(N) a model. Here we hope to clarify the 
situation. 

We use the notation of Ref. p| and consider (a 2 )^ = (0|a 2 |0) — iVm 4 , where a(x) is 
the auxiliary Lagrange- multiplier field with vacuum expectation value \fNrn 2 . (a 2 )^ is the 
next-to- leading-order correction xal/N and can be written as || 

d 2 k r——: r^r/, Vk 2 + Am 2 + ^¥ 



(a 2 ) « = -4tt / -—J- Jk 2 (k 2 + 4m 2 ) I In 

Jk<a ( 2TC ) ' 



rA{n) ( x - 1)2 

-m 4 / dx K - ! , 9 
Ji x z mx 



where 
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Am 2 2m J 



(10) 



We would like to point to the fact that the integrand in Eq. (|9]) has no infrared sin- 
gularities even though its large k behavior contains a factor of the running coupling 



f{k 2 ) = (47r)/ ln(A; 2 /m 2 ), which has an infrared pole at k 2 



m 



As discussed below, 



the QCD gluon condensate has a similar feature. Not surprisingly, (a 2 )^ shows a power 
divergence for large [i. Thus the integral in Eq. (|9|) has no way of producing any ambigu- 
ous imaginary part in contrast to Ref. f7j. A direct integration shows (a 2 )^ contains the 
exponential-integral function Ei(ln(/i 4 /m 4 )) for fi 3> m. Here one does not make any 
assumption or prescription in performing the integral and Ei(x) is a well-defined principle- 
value integral with a known asymptotic behavior 



Ei(x) = -- \Ex(-x + ie) + E x (-x 



~ e 



n=0 



nl 



X 



71+1 



{x > 1) 



(11) 



The only ambiguity that comes in is when all we know is the asymptotic series in Eq. (|TTJ) 
and we have to use a Borel transform to resum it, in which case the Borel transform has a 
singularity on the positive real axis; this is the infrared renormalon. Note also that, although 
the series in Eq. (O) is not unambiguously Borel summable due to the infrared renormalon, 
it is still a useful asymptotic expansion since its first [x] terms, with [x] meaning the integer 
part of x, is still a good approximation to Ei(x), as can be checked numerically. Thus if 
we allow power terms to include powers of /i multiplied by an asymptotic series (Borel- 
summable or non-Borel-summable) in the coupling constant, the separation of power terms 
from the nonperturbative contributions is possible. In the present case, the expansion of 
m 4 Ei(ln(// 4 /m 4 )) provides an example of such power terms. 

Now we apply these ideas to the case of the photon condensate 

(Hn)F 2 ) = (0\(a(^)MF^\0) , 

where a(/i 2 ) = e 2 /47r is the fine structure constant renormalized at /i and is the elec- 
tromagnetic field strength tensor. We note an earlier discussion on the photon condensate 
in Ref. plfl . We can define the photon condensate similar to Eq. ([Tp with <p(x) replaced by 
F^ v (x). Note that this is a gauge-invariant definition. With this definition, one can convince 
oneself that 



[-F 2 ) = -a( l x 2 )i l i A -° 



I (^ ^ ~ k ^°^ 



(12) 



where some scheme- dependent subtraction is understood, the momentum is Minkowskian 
with the metric g^ v = diag(l —1 —1 —1), and D liv {k) is the full photon propagator renor- 
malized at fi. In a general covariant gauge as represented by the parameter A, current 



conservation allows one to write |22 
D 



(A;) 



-9fw + 



k^kjj 



k 2 



1 



k^k u 



l + n(-fc 2 ) A k 2 



k 2 



(13) 



where U(—k 2 ) is the photon polarization. Plugging Eq. (|l"3"D into Eq. (|T^) and Wick-rotating 
to Euclidean space, we obtain a result independent of the gauge parameter, 
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\ ' tt^ J {2n) d l + n(F) ' 1 J 

Note that we cannot set the dimension of any integrals implicit in the expression for 
U(k 2 ) to four before we perform the explicit integration in Eq. fllij) . Keeping I1(A; 2 ) at 
one-loop order and evaluating Eq. ([[4]) order-by-order in a(fi 2 ) in DR, one finds 

<-F 2 > = . (15) 

7T 

However, it is instructive to perform the calculation by introducing an ultraviolet Euclidean 
cutoff (or subtraction point) \i. To maintain gauge invariance, we should use the result for 
U(k 2 ) from a gauge-invariant renormalization. Also power terms must be subtracted from 
the final result as discussed earlier to compare with the result in DR. 

Once we set d = 4, the integrand in Eq. fll4|) becomes the running coupling constant 
a(k 2 ). Also, the integral is dominated by the region of momentum k ~ /i, since the phase 
space provides a factor of k 2 , which suppresses the momentum region of very low momentum. 
For n the electron mass, we can then use the asymptotic form of the running coupling. For 
our purpose here, it suffices to keep only the one-loop part which corresponds to summing 
up the one-loop bubble diagrams. This is the well-known expression 

= «Q 2 ) 3tt 

^ ' l-K/i 2 )/3vr]ln(P/^) - ln(P/A2 ED ) ' 1 > 

where A 2 _ the location of the Landau ghost pole. After a Wick rotation, 

the integral in Eq. (|14T ) can be written in terms of the exponential-integral function. For 
/i <C Aq ED , it equals [^U 



9 rA 4 QED Ei(ln(A 4 QED ///)) = ^ £ , (17) 



g7r 2 QED 1 V*~V QED/ r ^3 ^ ^ 

which are all power contributions and must be subtracted to compare with the result of DR. 
Thus we arrive again at a vanishing photon condensate. 

We can deal similarly with the gluon condensate of the physical QCD vacuum |0): 

(^Gfl) = {0 \<^l G a G ^a l0) f a = l,2,-..,N 2 -l (18) 

7T 7T 

where a s (—k 2 ) = g 2 {— k 2 )/{Au) is the running strong coupling and the gluon field strength 
tensor is 

G% = d,A a v - d u Al + gf abc AlAl , (19) 

and 

G^G™ = 2{d il A° v - dvA*) d"A ua + Agf^d^A" A\A C V + g 2 f abc f ade A\A C V A^ ' A ve , (20) 

with f abc the structure constants of the SU(N C ) color group. In contrast to the photon 
condensate, a definition analogous to Eq. (JTJ) may not apply since T[G a ^ y (x)G^ va (ti)] is not 
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gauge invariant. Furthermore, mixing between the G 2 operator and gauge- variant operators 
must be taken into account if a gauge with ghosts is chosen [|4J]. Nevertheless, it is tempting 
to generalize the QED result in Eq. (|I4| ) to obtainF] 

71 71 J [Z7l) a 

where the (N 2 — 1) factor comes from the sum over colors. Note that Eq. fl2~T|) involves the 
running coupling at any given scale to all orders and the leading-order term coincides with 
the result in Ref. 0. 

To give some justification for Eq. (pT|), we can include the G 2 operator in the QCD 
lagrangian with a constant Lagrange multiplier £ and choose a physical gauge such as the 
axial gauge with gauge fixing term 

-±.(n.A a ) 2 , (22) 
where n M is some arbitrary fixed vector. Since the ghosts decouple |25|j26[] , there is no 



operator mixing with the gluon operator in the chiral limit. After evaluating the usual 
connected generating functional as a function of £, we can obtain the gluon condensate 
by taking a derivative at £ = 0. Note that it is convenient if one rescales the fields by 
letting A — > Aj y/1 + £, and similarly for g and A. This way one can show that, at least 
at one-loop order, the gluon condensate is independent of the gauge parameter because g 2 X 
is not affected by the above scaling. The result turns out to be the same as that obtained 
from using Eq. ([!]) with <p(x) replaced with G*(x) and gauge-parameter-dependent terms 
disregarded. As a result, the first term of Eq. (|20D alone will produce the result of Eq. ([H]) at 
least at one-loop order if one uses the known one-loop gluon polarization in physical gauges 



25| ) [27j| . Note that, in perturbation theory, one can show order-by-order in dimensional 
regularization that all the terms in Eq. (|20|) yield vanishing contributions to the condensate, 
although the nonperturbative sum to all orders may not vanish. Intuitively, we expect that 
the first term of Eq. (|20|) should dominate the nonperturbative gluon condensate since it is 
a one-body operator and its matrix element involves only the creation and absorption of a 
single gluon while the matrix elements of the rest of the terms involve two or more gluons. 
Thus we conjecture that Eq (ETJ) is a good approximation in QCD. 

It is interesting to see what we can learn from Eq. (pi] ) for the perturbative QCD vacuum, 
which is here defined by requiring the running coupling to be given by the asymptotic- 
freedom result. At the one-loop level, we have 



(K( '' J) = PM(k 2 -ie)/Al co ) ' (23) 



where (3 = (11N C — 2N{)/3 is the first coefficient of the /5-function of QCD with Nf light 
quarks and A QCD is the QCD scale parameter. Note that we have retained the Feynman ie 
prescription even after the Wick rotation to Euclidean space. 



2 We note some estimates of the gluon condensate from the gluon propagator in Ref. 
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We again use a cutoff regularization with ji > A QCD . In earlier work, integrals similar to 



Eqs. (21) and (23) have been examined to show the factorial growth of the coefficients of 



perturbation theory f28 |. Here we note that our specification of the Feynman ie prescribes 



the way to go around the infrared singularity and allows to perform the integral in Eq. (|2T|) 
analytically. The principal value of the integral gives a real part that can again be written 
as an exponential-integral function and integration around the infrared singularity leads to 
an imaginary part. We find 



where we have subtracted the power terms: 
3(iV c 2 -l) A4 TO , 1 , 4/A4 „ 



, k2 , A 4 QCD Ei(ln(/, 4 /A 4 QCD )) = j^{Nl - l)E;!©"W^t +1 ■ (25) 

Thus the energy density of the perturbative vacuum is imaginary, indicating its instability. 

One may ask: where are the ambiguities of the infrared renormalons? As discussed 
earlier, the renormalon ambiguities appear only when one tries to resum the asymptotic 
perturbation series in Eq. (p5|) using a Borel transform. The Feynman ie prescription gives 
the same result as that obtained from the Borel summation method when the integration 
contour in the Borel plane is deformed to pass above the infrared renormalon pole. Thus 
the Feynman ie prescription says unambiguously that the perturbative QCD vacuum is 
unstable as expected. Notice also that, had we considered the power terms in Eq. fl25|) as 
the perturbative contribution to the phenomenological gluon condensate, at fi ~ 1 GeV this 
contribution would dominate the phenomenological value of ~ (350 MeV) 4 determined from 
QCD sum rules. In Ref. Q, the estimated power contribution is small because a low scale 
\i is used. 

At present we know little about the physical running coupling constant of QCD in the 
low-momentum region, although we expect a s (k 2 ) to exhibit asymptotic freedom for large 
momentum. However, given the phenomenological value of the gluon condensate from QCD 
sum rules, we can consider Eq. ( |2T| ) as a constraint for the running coupling. Such a 
constraint may be important in deciding what Ansatz to use for the running coupling in 
some QCD phenomenology such as the DSE approach. 

In the literature, the Ansatze used in the DSE approach are a combination of the fol- 
lowing: 

a ° (k2) = /Wr + tf/AScJ ' (26) 
ai (k 2 ) = 7tak 4 e- k2/A , (27) 

= ^V fe2 / A , (28) 

where r, a, A, and x are phenomenological parameters. These yield the following corre- 
sponding contributions to the condensate: 
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;f G\ = ^l^ 1 ' Kou [r Ei(lnr) - Ei (lnr 2 )] , (29) 

.^ l= m^i^ t ( 30) 



4vr 2 



<- G2)2 .£M^> x * A . (31) 

7T l07r 2 

Notice that ao(& 2 ) is an extension of the one-loop renormalization group result to low- 
momentum region; typically r is taken to be 2 to 10. To obtain Eq. (ffil), we have sub- 
tracted power terms, which are the same as Eq. (|25|) for /i 2 ^> A 2 CD . For N c = Nf = 3 and 
r = 3, we find ((a s /7r)G 2 ) = (0.568A QCD ) 4 , while, for r = 10, ((a s /7r)G 2 ) = (1.44A pCD ) 4 . In 



Ref. the running coupling is taken to be a s (k 2 ) = ceo(k 2 )+ai(k 2 ) with a = (387 MeV) 4 , 
A = (510 MeV) 2 , r = 10, and A QCD = 228 MeV. The resulting gluon condensate is 
(788 MeV) 4 . In Ref . the running coupling is a s (k 2 ) = a (k 2 )+a 2 (k 2 ) with x = 1.14 GeV, 



A = 0.002 GeV , r = 3, and A QCD = 190 MeV. This gives a gluon condensate of (191 MeV)' 
Thus both results differ significantly from the value determined from QCD sum rules. Nev- 



ertheless, the parameters are somewhat flexible, as indicated by the authors in Ref. p3 
Thus our result provides a further constraint on the parameters. 

Finally, we remark that the gluon condensate as given by Eq. (BT) is not sensitive to the 
behavior of a s (k 2 ) near the origin because of phase-space suppression. Thus, if confinement 
is associated with the running coupling near the origin, our result indicates that the nonzero 
gluon condensate is not neccessarily related to confinement (see also Ref. ||23|| ). 



To reiterate, we have argued that the power terms can be excluded consistently from the 
OPE. We have shown that the photon condensate vanishes and suggested that the gluon 
condensate can be simply related to the physical running coupling of QCD. Finally, we 
examine some Ansatze for the running coupling that have appeared in the literature. The 
resulting values for the gluon condensate are different from the value determined from QCD 
sum rules. 
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